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ABSTRACT 


Digitized cloud brightness data for July 1969 are examined to 
determine the degree of relationship between organized brightness 
patterns and the large-scale convection field in the tropical 
western North Pacific. Correlation coefficients between kinemati- 
cally computed vertical motion fields and cloud brightness are 
generally low, except for a few notable days. This may be due to 
the quality of the available data, particularly the vertical motion 
field. Nevertheless, indications of better correspondence between 
the two fields are noted late in the month for the western portion 
of the region, for the latitude band between 10°N and 20°N. 

Examination of time-longitude sections reveals a close associa- 
tion between propagating brightness patterns and vertical motion 
fields, indicating that the former are reflections of synoptic wave 
passages, rather than simply inactive clouds advected by zonal flow. 

In view of the potential usefulness of such satellite data, a 
technique iS proposed that uses satellite cloud data to objectively 


determine the large-scale tropical flow. 
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L. deeepuCTIOCH 


Tropical analysis has long suffered from a sparse observational 
network. It is natural to hope to supplement this network with avail- 
able satellite data. Digitized cloud brightness data, acquired by the 
Advanced Vidicon Camera Systems of the operational satellite system, 
have been available since 1967. Methods of utilizing this data on an 
operational basis may contribute significantly toward improving 
synoptic-scale analysis of the tropics. 

Satellite photographs of tropical regions reveal that much of the 
tropical cloudiness is in the form of large, connected masses of very 
bright clouds. These masses are composed of many individual cumulus 
clouds, whose extensive cirrus canopies account for their brightness. 
Through the use of satellite photographs, these large, bright cloud 
masses, or “cloud clusters" as they are now called, have received con- 
Siderable attention in recent years. 

Chang (1970), using time-longitude series of satellite photographs, 
found that many cloud clusters over the tropical Pacific are well 
— and have characteristics resembling those of synoptic-scale 
westward propagating wave disturbances. Reed and Recker (1971) showed 
further evidence that the propagating cloud clusters are linked to 
synoptic-scale waves in the wind field, while Williams (1970), using a 
compositing technique, examined the mean structure of trade wind cloud 
clusters in the western North Pacific. Wallace (1971) also found a 
high degree of organization to cloud brightness patterns, with many 


patterns maintaining their identity over a period of weeks. 





Work such as this has indicated that clusters represent the or- 
ganizea cumulus convection whose importance to large-scale tropical 
motions is well recognized. On the belief that clusters do represent 
organized convection, attempts to utilize cloud brightness, either 
directly from satellite photographs or in a digitized form, have 
Mainly taken two directions: 1) to estimate rainfall, and 2) to 
indicate vertical motion. 

Precipitation, while important in its own sake, also indicates 
the spatial extent of convection and updraft intensity, as well as 
latent heat release. Attempts to estimate rainfall using satellite 
data have been reviewed by Martin and Scherer (1973). Specifically, 
Woodley and Sancho (1971) and Woodley et al (1972) have attempted to 
relate rainfall to cloud brightness using sophisticated methods to 
correct or normalize the raw cloud brightness data before using it 
quantitatively. . 

In examining digitized cloud brightness data, Wallace (1971) 
found a tendency for it to be well correlated with certain indicators 
of disturbed weather, particularly vertical motion, and he emphasized 
his belief that digitized brightness data may be the best availabie 
paiatiox ofvdisturbed weather dm the topics: Wallilace and Ciatig 
(1972) applied digitized cloud brightness data to the study of tropi- 
cal wave disturbances. They found fairly good correspondence between 
area-averaged cloud brightness and vertical motion, and expressed 
their belief that it may be possible to obtain = estimate of the 
vertical motion field over the tropics from the cloud brightness data 


alone. 





If, as these studies suggest, a vertical motion field over the 
tropics can be obtained from brightness Gata alone, then aemmayebe 
possible to derive the large-scale tropical flow field even over 
regions where wind observations are extremely Sparse. 

This study is an attempt to empirically relate the organized 
brightness to the large-scale convection field through the vertical 
motion field. The objectives of this study are 1) to examine the 
relationship between digitized cloud brightness and the large-scale 
convection field in the tropics, and 2) to propose a technique to 
improve the analysis of data-void tropical regions through the use 
BP scaigitized cloud brightness data. 

The area selected for this study is the tropical western North 
Parelnic, Lrom the equator to SO°N, I25°E to 180°. =the regicen an 
cludes the Caroline and Marshali Islands with nine rawinsonde sta- 
tions. Figure 1 shows the region and the rawinsonde network. 

Digitized cloud brightness and upper-level wind data used in 
this study were provided by the National Center for Atmospheric Re- 
search (NCAR). The cloud brightness data, acquired by NCAR from 
the National Environmental Satellite Service (NESS), are digitized 
meiies averaged for 5° x 5° latitude-longitude squares. The units 
of brightness are arbitrary, but have been designed to show a 
linear rate of change from black to white within a scale of zero to 
ten. (This data, although an average, will be referred to as raw 
cloud brightness data in this paper.) Upper-level component winds 
at the 700, 500; 300, 250 and 200 nib levels were acquared by Nek 


from the National Meteorological Center (NMC) tropical grid analyses. 





These winds are also on 5° x 5° latitude-longitude grids, and are the 
result of analyses based on all tropical rawinsondes, wind soundings, 
aircraft observations, and some few winds deduced from cloud drift as 
seen from satellite photographs. NCAR also provided rawinsonde obser- 
Vations for the region of study, which were used to obtain surface 
data and low-level wind fields. 

The time period of the study is July 1969, which contained 23 days 
of useful data. The O000 GMT wind data were selected as being closest 
to the time of the satellite observation. The satellite in operation 


Sermo July 1969 was ESSA Ix. 
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Ti RELATIONSHIP SBETWEENMGLOUD BRIGRINESs 
AND THE LARGE-SCALE CONVECTION FIELD 

To determine the relationship between cloud brightness and the 
large-scale convection field, the digitized cloud brightness data 
is compared both subjectively and objectively to the vertical p-veloc- 
ity field (Ww). The vertical motion field 1s chosen because Wallace's 
(1971) results indicate the potential for good correlation. The raw 
cloud brightness data as available from NESS is used with no attempt 
to adjust it for radiation geometry considerations or instrument 
characteristics. 

Fields of vertical motion W and divergence are computed kinemati- 
cally. Divergence computations are performed at the 9300, 700, 500, 
300, and 200 mb levels by simple finite-differences over the 5° x 5° 
grid. Vertical velocities are constrained by the boundary conditions 
W = O at 1000 and 100 mb. The correction to satisfy these boundary 
conditions on W is distributed linearly with pressure to result in 
adjusted divergence values at the stated levels, and values of WwW at 
the 900, 800, 600, 400, 300, and 200 mb levels, in addition to the 
boundary values. The validity of such calculations ee a 5° & S2oreid 
1S undoubtedly questionable. However, for the purpose of correlating 
the vertical motion with cloud brightness, an approximate indication 
of the vertical motion should be adequate to at least show whether the 
correlations tend to be high or low. Furthermore, the computed 
divergence profile in areas of organized brightness shows a deep layer 


of convergence with a thin layer of divergence above, and both the 
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divergence profiles and the vertical motion profiles in areas of or- 
ganized brightness show generally good agreement with the results found 
by Williams (1970) and others, thereby lending some confidence to the 
calculations. Nevertheless, in view of the calculations, more confi- 


dence is placed in the directly observed fields of cloud brightness. 


As SMOOTHING OF CLOUD BRIGHTNESS FIELDS 

Intuitively, the best correspondence between cloud brightness and 
vertical motion should occur when the clouds are well organized, active 
and deep -~ indicating strong rising motion through much of the tropo- 
sphere. The results of Woodley ct al (1972) support this. In examining 
the correspondence between radar echoes (which indicate precipitation 
and therefore probable areas of deep convection) and cloud brightness, 
they found the best poe cenoneence when clouds are young and active and 
the worst when clouds are old and dying. The brightness patterns should, 
therefore, be maintained in both time and space. Small, isolated 
bright spots will be undesirable, as well as dissipating brightness 
patterns no longer supported by convection. To separate small, isolated 
bright spots from the organized convective brightness, hopefully retain-~ 
ing only brightness patterns representing cloud clusters, the brightness 
data . subjected to several elementary smoothing techniques. It is 
found that simple area-averaging of surrounding grid values with various 
weighting applied to the central value is adequate to emphasize the 


larger, Organized brightness areas and wipe out isolated bright spots. 


BY. SUBJECTIVE COMPARISON OF CLOUD BRIGHTNESS TO VERBICAL MOTIONS FSELD 
The unmodified and smoothed brightness fields are analyzed, or con- 


toured and then compared with the analyzed vertical motion fields at 
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the 600 and 300 mb levels, which were selected as representative of 
the low and upper-level fields, respectively. The brightness fields 
show good organization, with many brightness patterns maintaining their 
identity as they move westward across the region of study. The verti- 
cal motion fields show alternating patterns of rising and sinking 
motion that propagate westward witn time and have spatial characteris- 
tics typical of synoptic-scale wave disturbances. This apparent wave 
structure becomes most evident at the 300 mb level during the latter 
part of the month. These wave disturbances and their relationship 
to the organized cloud patterns will be discussed in more detail later. 
The best correspondence between brightness and vertical motion 
appears when the rising motion is deep, extending up through the 300 
mb level. Accordingly, in most cases, the 600 mb comparisons reveal 
nothing that is not apparent at the 300 mb level. In fact, the rela- 
tive orientation and size of the patterns appear better at the 300 
mb level. 
Figures 2, 3 and 4 show the contoured fields of 600 mb vertical 
motion, 300 mb vertical motion and cloud brightness, respectively, 
mor 2/ July 1969. Figures 5 and 6 show the 300 mb vertical motion and 
the cloud brightness Eield, respectively, for 29 July 1969. The 
figures for 27 July compare well, and in fact, exhibit the best corres- 
pondence encountered during the month. The figures for 29 July are 
more typical. Note that there are essentially three areas of strong, 
organized brightness, none of which compare well to the vertical motion 
field. The contoured cloud brightness fields shown are fields of raw 


brightness data. While the smoothed brightness fields are more pleasing 
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to the eye, the raw brightness fields are found to lend themselves 
more readily to visual comparison with the wertical moeron frelds- 

Comparison of contoured fields for the entire month reveals the 
following: 

a. Better correspondence between areas of organized brightness 
and centers of rising motion is seen for the western portion of the 
region and between the latitudes of 10°N and 20°N. In fact, the 
eastern portion possesses comparatively few organized brightness 
patterns, exhibiting generally a widespread area of zero or low 
brightness. ThiS increase in organization and intensity may be a 
reflection of the westward intensification of cloud systems found 
by Chang et al (1970) and Reed and Recker (1971) as wave disturbances 
move across the region of study. Sikdar et al (1972) also noted 
marked geographical variations in cloud activity over the Pacific. 

b. The correspondence tends to improve toward the end of the 
month. Thus, the organization of the brightness patterns may vary 
considerably even within a season, or there may be a steady increase 
in the organization through the late summer. 

c. While Wallace's (1971) results indicated a good phase relation- 
tks between brightness and vertical motion, this is not evident in 
this study. However, the brightness centers are displaced both east 
and west of the centers of rising motion, which suggests that the fault 
may lie in the inaccuracies of the vertical motion calculations. Fur- 
thermore, there is a time difference between the two fields which is 
ignored in this study. The vertical motion field is computed from 


OOOO" GMT data. The time of the Satellite CbServalion Varies across 
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the region, but generally is within two to three hours of the 0000 
GMT observations. 

a. Many of the brightness patterns exhibit a southwest to north- 
east tilt, which seems to correspond to the usually observed tilt of 
synoptic waves in the tropics, as summarized by Wallace (1971). 

e. Except on certain days and for limited regions, the corres- 
pondence between areas of brightness and areas of rising motion cannot 
be considered satisfactory. In general, for this study, centers of 
organized brightness may appear in regions of either rising or sinking 


motion. 


Seo reneCTIVE COMPARISON OF CLOUD BRIGHTNESS TO VERTICAL MOTION FLELD 

To objectively determine the relationship between cloud brightness 
and the vertical motion field, scatter diagrams and simple linear cor- 
relation coefficients are obtained. Due to the poor quality of the 
computed vertical motion fields, the objective is not necessarily to 
seek high correlation coefficients, but rather to determine if the 
correlation tends to be generally high or low. 

Vertical motion fields at the 800, 600, 400, and 300 mb levels are 
used for the correlations with both the raw and smoothed brightness 
fields. The eee re are computed on both a daily basis and for 
the month. In addition, on the basis of the subjective comparisons, 
correlations are obtained between the vertical motion and the smoothed 
brightness for just the western portion (125°E to 155°E) of the region, 
and for only the last ten days of the month. 

Since the subjective comparisons suggested that the best corres- 


pondence exists when the rising motion is deep, the correlations are 
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expected to be generally negative at all levels. That iS, increaSing 
cloud brightness as the vertical motion w becomes more negative (stronger 
rising motion). This 1S indeed the case. The highest correlation is 
obtained most frequently at the 300 mb level. In fact, in 92 percent 

of the cases, the best correlations occur at the upper levels, at either 
BUO-or 400 mb. 

Figures 7 and 8 show scatter diagrams of smoothed cloud brightness 
versus 600 and 300 mb vertical motion, respectively, for 27 July 1969. 
Figure 9 shows the scatter diagram for the raw cloud brightness versus 
the 300 mb vertical motion on 27 July 1969. As mentioned earlier, this 
day illustrates the best correspondence found during the month. The 
correspondence is best at the 300 mb level, and for the smoothed bright- 
ness case. While the scatter 1S greater at 600 mb, the correspondence 
is still reasonably good. The linearity is particularly good for the 
smoothed brightness data. 

A more typical case is shown in Figures 10 and 11, which are scatter 
diagrams of smoothed cloud brightness versus 600 and 300 mb vertical 
Mouron, respectively, fou 29 duly 2269.. Note a great dealmere Geatter 
at the 600 mb level than. at the 300 mb level, and note also that the 
Per at the 300 mb level 1s about a line of nearly zero slope. In 
other words, the entire range of brightness values exists with both ris- 
ing and sinking motions. This situation occurred quite frequently. 

Another common situation is illustrated by Figure 12, wnich shows 
the scatter diagram of 300 mb vertical motion versus smoothed cloud 
brightness for 22 July 1969. There appear to be two groups of bright- 


ness values, one group of low brightness values and another of high 
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values with the split at around a value of 2.0. This suggested that 
establishing a "brightness threshold" by deleting brightness values 
less than a mean brightness value or an arbitrary value around 2.0, 
from the cloud brightness sample, might improve the correlations. 
However, aS will be discussed later, such modifications do not actually 
improve the correlations. 

Table I lists the correlation coefficients between the 300 mb ver- 
tical motion and the raw cloud brightness, the smoothed brightness, 
and the smoothed brightness for just the western portion of the region. 
This level produced the highest correlations. The coefficients are 
listed by day, for the month, and for the last ten days of the month. 
The following comments apply: 

a. For the full grid, the daily sample size is 72. Panofsky and 
Brier (1965) describe an analysis of variance method to determine if a 
€em@eelation coefficient is significant, i.e., if it is greater than 
eneeewhich could occur by chance alone. Using thas method, the corre= 
lation coefficients must be greater than about 0.24 to ensure signi- 
hieence at the 5 percent level and 0.30 at the 1 percent level. On 
mars basis, ten of ae 23 Gays for the smoothed brightness case show 
significant GG@rcelations abkithe oS »percent level, 2nd.fivercats tie sla pen— 
cent level. me the western portion of the grid the daily sample size 
1s 42. For this sample, the coefficients must be greater than 0.30 to 
ensure significance at the 5 percent level, and 0.38 at the 1 percent 
Mevel. Thus, for the western portion, 11 dayswshow significant co- 
efficients at the 5 percent level and eight days at the 1 percent level. 
Although this difference is certainly small, it does suggest (as did 


the subjective comparisons) that generally the western portion of the 
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region for the area-averaged brvgimedess has) a higher dequcevor ia. 
tionship to the vertical motion field. 

b. The correlations are negative on all except two days, the first 
and the twelfth. These are days on which there are few centers of 
EiEsing motion. 

c. Subjectively, it is easier to compare brightness with vertical 
motion using the raw cloud brightness fields. Objectively, while 
smoothing in most cases improved the correlations, the improvement is 
quite small. 

dad. As suggested by the subjective comparisons, the correlations are 
higher during the latter portion of the month. (For the western sample, 
nine of the last ten days have coefficients which are significant at 
the 5 percent level of confidence.) 

Thus, generally the correlation between organized brightness patterns 
and the vertical motion field for this study is guite low (only about 
half of the days show coefficients higher than that which could occur 
by chance alone). However, as noted above, there is a tendency for 
higher correlations te occuUmMGor certain réegiens and time peéeiods. 

Correlation coefficients were also computed for several modifica- 
tions. For example, a "brightness threshold" was established. Bright- 
ness values less than the monthly mean of cloud brightness and 
arbitrary brightness values of one and two were deleted from the 
smoothed brightness sample. The remaining samples of brightness values 
were then correlated with the vertical motion Field, and also with only 
the field of rising motion. In all cases, such modifications resulted 
in only slightly higher coefficients but reduced the sample size to one 


of questionable significance. 
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D. COMPARISON OF Titis-LONGITUDE SECTIONS 

As mentioned earlier in this Section, the analyzed fields of ver- 
tical motion exhibit wave disturbances which become particularly 
evident at the 300 mb level during the latter portion of the period 
of study, and many of the brightness patterns exhibit a tilt in the 
same orientation as that of synoptic-scale waves observed in the wind 
field. These results suggest that both fields may be organized by 
synoptic-scale waves. Thus, one may expect to find a better relation-~ 
ship between the two fields in a time-space domain corresponding to 
that of synoptic~scale waves. If this is the case, it may also pro- 
vide an answer to a question addressed by Chang (1970) and Wallace 
(1972) regarding the westward propagating cloud patterns. 

Chang (1970) and Reed and Recker (1971) suggested that propagating 
cloud patterns may be reflections of westward propagating synoptic-~ 
scale waves. However, Wallace (1972) noted that the westward propa- 
gating cloud patterns seen on time-~composite satellite photographs 
generally move with nearly constant speed, but that the horizontal 
spacing between such clouds appears to be quite variable. In view of 
this "non-dispersive" phase speed, he expressed a doubt that these 
Beads are reflections of wave passages, but that, more likely, they 
are simply advected by zonal flow. 

In order to examine these questions, time-longitude sections of 
the raw brightness and the fields of rising motion are compared for 
five-degree latitude bands. Figures 13 and 14 show, respectively, 
the sections of the cloud brightness and the rising motion for the 


lo = 15°N latitude band. This latitude band isthe most userul tor 
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comparisons as other bands show fewer propagating cloud patterns. The 
path of westward moving brightness patterns can be quite easily traced 
from these sections. Several such paths are plotted as propagation 
lines in the brightness section (Figure 13), and then overlapped on the 
vertical motion section (Figure 14). It can be seen from the latter 
section that the cloud pattern generally propagates along an area of 
strong upward motion. Considering the inherent inaccuracies in the 
two fields and the time differences, this result may be viewed as a 
strong indication that the propagating cloud patterns are indeed asso- 
ciated with synoptic wave passages. This is due to the Pace that sre 
they are simply inactive clouds advected by zonal flow, one would not 


expect such an associated upward motion tendency. 
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III. DESIGN OF OBJECTIVE TECHNIQUE THAT 
USES DIGITIZED BRIGHTNESS DATA 

Although a firm relationship between satellite-observed cloud 
brightness and meteorological parameters has yet to be established, 
confidence must be placed in what is actually seen through the satel- 
lite camera. The viSual indications that the dense cloud masses 
reflected in satellite photographs are related to active convection 
and large condensation heating can hardly be ignored. With the 
belief that such a relationship exists, and in view of the need and 
desirability of using satellite data to improve the analysis of a 
data-void tropical region, an objective technique has been designed. 
Although the model is still in the development stage, undergoing 
further modifications and testing, the basic design is outlined here 


to illustrate such a possible technique. 


A. BASIC MODEL 

The baSic idea of the model is to consider that the horizontal 
divergence term is the forcing function in the vorticity equation, 
and ie assume that the unknown divergence field may be replaced by a 
satellite-observed brightness field. The linearized vorticity equa- 
tion is then integrated to yield a steady state perturbation stream 
function resulting from such forcing. 


The vorticity equation in pressure coordinates may be written as 


3 3 ; a dw 
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where 
t = the relative. vorticity 
u = zonal component of the non-divergent flow 
Vv = meridional component of the non-divergent flow 
f = Coriolis parameter 
8 = df/dy 
dW/dp = vertical differentiation of the vertical p-velocity 
V = total horizontal wind vector 
V = isobaric gradient operator 


The thermodynamic energy equation is 


where T is temperature, W the vertical p-velocity, S the static sta- 
bility, and 0 the diabatic heating rate. Due to the usually small 
temperature fluctuations in the tropics, the first two terms of 
equation (2) may be neglected in comparison to the heating term. 
Thus, 


WS = QO (3) 


which is a statement of balance between diabatic heating and adiabatic 
Cooling. Stbstituting (3) into the right hand side of (1), the 


divergence at a level becomes 
a | (4) 
where oF is a measure of the magnitude of the heating, and dn/dp repre- 


sents the vertical gradient of the heating profile. The parameter oh 


Ze 





is assumed to be proportional to the satellite-observed cloud bright- 
ness and thus may be empirically related to the brightness field. The 
term on/op, for simplicity, is specified as a constant based on obser- 
vational studies by Reed and Recker (1971), Williams (1970), and Wallace 
(1971). Expression (4) may then be viewed as the known forcing for 

ene vorticity equation (Chang, 1973). 

For initial testing of the model, the simplest and most direct 
method is to relate the brightness empirically to the 300 mb divergence 
fHreld by a proportionality factor. This brightness forcing may then 
be substituted for the horizontal divergence term in Equation (1). 

Assuming that the basic state of each parameter is a function only 
of time and the y-axis, the linearized vorticity equation may be 


written as, 


2a 


fee, - oO 2. OW’ Os, OW ar 
ave Vew'+ u Vw a +B = Gr ay) Bilton 13 — alee (5) 


ox dy? OX 
where primes denote perturbation quantities, bars denote basic state 
quantities, W is the stream function, B represents the brightness as 
related to Bier aencen and Dis a damping term. The damping term is 
included on the basis of a study by sHolton vandyGoltonw(19/2)scofsthe 
Meriicity panenes at 200 mb in the tropics, wherein they found that 
rather strong damping is necessary to achieve realistic stream function 
fields. This damping may be interpreted as due to the vertical trans- 
Bort of vorticity by cumulus convection. An sadditional term, ¢'B’, 
not shown in (5) may be included with no difficulty since it is still 


linear in computation if B' is specified. 
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The corresponding finite-difference equation is as follows, 





t-1 2 
(Vow. .) ~*{1+At (B-D)] a 
(V7. wie = ot ra a *{ ( aaa an 
td [b= Ae eee) | [1-At (B-D) ] 
72 _y2 _m _ ne a eee 
ee a an, 9g a hia, 5) a ee) 
+(- Bey -W ee, (ae er (6) 
a. 2m eed Epey fleet bre me 20) aaa Pool 


where the primes have been dropped for convenience, superscripts denote 
the time step, mis the map factor, d the grid distance, and the other 
terms are aS previously defined. This equation requires a mean zonal 
wind, whose determination is a major problem. However, it iS expected 
that monthly averages of observed data will prove adequate. 

Equation (6) may now be integrated over a data-void region to obtain 
a steady state stream function field for the region if appropriate 
boundary conditions for the stream function are given. Thus, applying 
the model to a data-void region requires a judicious selection of the 
region itself. In most cases, data-void regions are bounded (although, 
not ideally) by areas with adequate wind information, typically island 
groups or edges of continents. A data-void region may, therefore, be 
Pee ened as a rectangular area or “box" surrounded on all four sides by 
Beees Of arbitrary Size within which there is wind information. If 
stream functions on the boundaries of these data-rich boxes can be ob- 
tained, the values on the sides of these boxes which are adjacent to 
the data-void region may then be used as boundary conditions for the 
data-void region. 

To obtain stream functions on the boundaries of the boxes, which 


are adjacent to the data-void region, thereby obtaining boundary 
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conditions for the data-void region, a method suggested by Hawkins and 
Rosenthal (1965) to compute stream functions uSing available wind data 
is employed. In addition to the outer boundary of each box, the method 
also uses the immediate inner boundary. To avoid reducing the area 

of the box, and to increase the finite-difference accuracy, all fieids 
are first reduced from a 5° x 5° latitude-longitude grid to a 2.5° x 
2.5° latitude-longitude grid by simple interpolations. The Helmholtz 


equation 


> as 
V=k ° W+ Vx , (7) 


is used to separate the horizontal wind field into a purely rotational 
part (, the stream function) and a purely divergent part (x, the 


velocity potential). Equation (7) can be written in tangential coordi- 


nates aS 

oy 9X 

ae Sia mae: So 
and OW Z Ne 

ar ee Os fey) 


where n denotes the outward normal direction, and s the tangential 
direction. With these equations, the procedure to obtain boundary 
conditions is outlined in the following six steps: 
(i) Relax V*y =)“ v, using Xx = 0 on the outer boundary, to 
obtain Xx over the grid. 
(ii) Calculate dx/dn on the inner boundary using center~differences. 
(iii) Integrate Equation (8) trapezoidally around the inner 
boundary to obtain W along the inner boundary. 


(iv) Calculate dy/ds on the inner boundary using center-differences. 
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(v) Integrate Equation (9) trapezoidally from the inner to the 
outer boundary to obtain on the outer boundary. Corner 
values are filled in by interpolation from the surrounding 
Grideipeasats- 

(vi) Discard W on the inner boundary. 

To summarize the application of the model, a data-void region is 
specified such that it 1S a rectangular area or box, and is surrounded 
by four adjacent boxes which contain wind information. The Hawkins 
and Rosenthal procedure is then appliea to each of the adjacent boxes 
to obtain stream functions on the boundaries of these boxes. Stream 
function values on the sides adjacent to the data-void box, then become 
boundary conditions for the data-void region. The vorticity equation 
(Equation 6) is then integrated over the data-void region (using bright- 
ness data to represent the forcing) to yield steady~state stream func- 


tion fields for the region. 
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IV. SUMMARY AND CONCLUSIONS 


Digitized cloud brightness data for July 1969 Bee examined by sub- 
jective and objective means to determine the degree of relationship 
between organized brightness patterns and the large-scale convection 
field an the sssepical westemn NortheRacifiac. Dadidyeventical awebion 
fields for the region are computed kinematically uSing analyzed wind 
fields based on OOOO GMT wind obServations. Comparisons are then made 
between the satellite-observed cloud brightness fields and the fields 
Gi vertical motions 

The correlation coefficients between the two ficlds are low, in 
most cases. However, consistent tendencies of better correspondence 
between the two fields are noted: 

a. The correspondence iS almost always better at the 300 mb level. 
This is consistent with the assumption that brightness is related to 
deep convection. 

b. The correlation shows an improvement toward the latter portion 
of the month. In the- spectral studies summarized by Wallace (1971), 
late summer-early autumn Seems to be the most active season for tropical 
disturbances. Thus, the period selected for this study may be too early 
in the season. Further study extending to later months may yield hetter 
results. 

c. The correlation appears better in the westernmost portion of 
the region. Chang et al (1970) and Reed and Recker (1971) have shown 


the marked increase of wave amplitude as disturbances move westward 
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from the vicinity of the Marshall Islands. In general, the eastern 
portion of the region of study simply shows a lack of high brightness. 

Particularly in the latter part of the month, the computed verti- 
cal motion fields possess wave-like patterns. A time-lcngitude 
section study of the two fields indicates a general correspondence 
between propagating brightness patterns and centers of rising motion. 
This result suggests that the westward propagating cloud patterns 
noted in satellite photographs are reflections of wave passages, 
rather than advected cloud systems. 

Based on the assumption that a dense cloud mass is related to 
large latent heating, a vorticity equation model which may be used 
to diagnose large-scale tropical flow patterns using satellite cloud 
data was eine nce This model, which assumes a correspondence between 
the cloud data and the divergence field, is derived from simple 
thermodynamic arguments. 

Concerning the correlation between brightness and convection para- 
meters, major problems exist. First, the "ground truth" - the vertical 
motion field used in thes Sstudyois dit ficulitetovobtain vaccumate lyainom 
the available data. The kinematic method posses large errors, while 
the vorticity equation method, which 1S actually the basis for the 
objective scheme proposed, requires knowledge of the damping due to 
cumulus transport. (It may be possible to obtain this knowledge by 
experimenting with the objective technique described in Section III.) 
The thermodynamic energy equation approach to determine W must be 


ruled out because of the near balance between adiabatic cooling and 


diabatic heating. 
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Secondly, Woodley and Sancho (1971), among others, have pointed 
out that large, bright cirrus canopies resulting from cumulus outflow 
may remain for relatively long periods after the active convection 
has ceased. These bright but dead clouds no longer represent the 
active convection, but are not removed by the area-averaging process. 
Smoothing in the time domain iS necessary to remove these inactive 
clouds, and a successful objective technique may markedly improve the 
correlations. 

It was observed in this study that, while the centers of organized 
brightness may agree in position with centers of rising motion, the 
brightness perimeter may extend further over areas of sinking motion 
(low level divergence) than can be justified by errors in the vertical 
motion calculations. Conceivably, edges of organized cloud patterns 
may be advected or dispersed away from the area of rising motion to 
such an extent that the correlations are significantly lowered. 

As mentioned earlier, the vertical motion field and the brightness 
field were not observed at the same time. Although the two fields 
for the region of-.study are within two to three hours of each other, 
the extent to which this time difference, uncorrected, effects the 
correlations is difficult to ascertain. However, Martin and Suomi 
(1971) found that uncorrected satellite cloud photographs could be 
used to study time changes in convection within three to four hours 
of local noon. Thus, it may be that cloud photographs taken close to 
local noon (which is the case for the region of study) do sufficiently 


indicate large-scale convection. 
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The digitized brightness data as available from NESS is subject 
EO adjustments. Woodley and Sancho (1971) and Woodley et al (19/72) 
in attempting to use brightness to estimate rainfall have found that 
in order to use the brightness quantitatively, they must first cor- 
rect it--or normalize it, taking into consideration such things as 
cloud-satellite-sun geometry, satellite characteristics, and even 
printing techniques. These corrections require quite sophisticated 
techniques and are not, at present, being applied to the brightness 
data by NESS. 

Finally, this study did not utilize infrared data. It is antici- 
pated that future inclusion of infrared data should be useful in 
isolating the bright, cold-topped clouds which truly represent the 


organized convection. 
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TABLE I 


Correlation Coefficients of Cloud Brightness vs 
300 mb Vertical Motion 


Last ten days 


é IL 


Are: 


Smoothed Smoothed 
Raw Brightness Brightness Brightness 
Date PulleGrra Full Grid L2ee t= 155 
i suly 1969 Males, pe es eae 
Z = 05 -.17 708 
3 4 -.24 -.26 
6 = -.22 -.39 
8 -.12 -.15 -.10 
9 -.13 -.16 -.06 
10 -.10 -.14 -.08 
J -.08 -.09 -.01 
Iz 4 8) a2 aval 
13 =-.05 202 -.23 
is ~.08 =. 25 -.25 
iF/7 -.23 27 =222 
ing) -.23 -.23 -.32 
20 = 224 -.24 -.42 
21 =i 22 -.19 -.33 
22 -.30 -.33 -.30 
23 -.32 -.4] -~.45 
24 mee -.26 -.37 
25 = SZ -.30 =o Sis. 
26 -.36 -.40 ~,55 
27 eo) =O2 -.63 
28 -.34 -.33 -.51 
20 relies, -.30 -~.39 
Month ae = Ae =. 20 
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Time-longitude section of upward vertical motion Woycday) fer 10°%> 15°N 


latitude band, July 1969. 


Figure 14. 
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